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Abstract 

We compute the next-to-leading order QCD corrections to the off-diagonal ele- 
ments of the decay-width matrix T entering the neutral i?-meson oscillations. From 
this calculation the width differences AT and the CP violation parameters (q/p) of 
Bd and Bg mesons are estimated, including the complete 0{as) QCD corrections 
and the l/rrih contributions. For the width difference AP^ we agree with previous 
results. By using the lattice determinations of the relevant hadronic matrix ele- 
ments we obtain the theoretical predictions APf^/P^ = (2.42 it 0.59) x 10"'^ and 
AP^/Ps = (7.4 ±2.4) X 10~^. For the CP violation parameters, we find \ {q/p)d\ - 1 = 
(2.96 ± 0.67) X 10^^ and \{q/p)s\ - 1 = -(1.28 ± 0.27) x 10"^ These predictions are 
compatible with the experimental measurements which, however, suffer at present 
from large uncertainties. 



1 Introduction 



B physics plays an important role to test and improve our understanding of the Standard 
Model (SM). By exploiting the production of mesons with a large boost, experiments 
such as BaBar and Belle can provide accurate measurements of the decay time distribu- 
tions and hence valuable information for the Bd-meson lifetimes, CP violation and mixing 
parameters. The production of Bg, on the other hand, is out of the reach of the S-factories 
and improved measurements of the related physical observables will come from the Run II 
at Tcvatron and from the LHC. Recently, the measurement of the mass difference Arud, 
which controls the frequency of B^ oscillations, has been further improved. The world 
average is now ArUd = (0.502 ± 0.006) ps~^ [1]. For Am^ and for the width differences AF^ 
and AFs, instead, only weak limits exist at present. Theoretically, the width differences are 
suppressed by a factor ml /ml with respect to the corresponding mass differences. In addi- 
tion, AFs is predicted to be larger than AF^, the latter being doubly Cabibbo-suppressed. 
The present world average for the width difference of Bg mesons is [1] 

AF,/F, = 0.07 , (1) 

where the constraint r^Bg) = T{Bd) has been used. Concerning the B^ system, a pre- 
liminary result for AF^/ has been recently presented by the BaBar collaboration [2]. The 
experimental result reads ^ 

AFd/Fd = 0.008 ± 0.037(stat.) ± 0.018(syst.) . (2) 

In ref. [2], a preliminary measurement of the parameter \{qlp)d\i which defines the 
amount of CP violation in the mixing, is also presented 

\{(l/p)d\ - 1 = 0.029 ± 0.013(stat.) ± 0.011(syst.). (3) 

In the SM, the CP violation parameters \{q/p)d\ — 1 and \{q/p)s\ — 1 are suppressed, with 
respect to the width differences, by a factor m^/ml. Moreover, \{q/p)s\ — l has an additional 
suppression factor (where A is the sine of the Cabibbo angle) with respect to \ {q/p)d\ — 1- 

The phenomenology of Bg — Bq oscillations is described in terms of a 2 x 2 effective 
Hamiltonian matrix, M'^ — iV^ /2. The b-quark mass, being large compared to Aqcd, allows 
to apply an operator product expansion (OPE) to the calculation of the decay- width matrix 
F^ [4, 5] . In the case of the off-diagonal matrix elements, which the observables AF and 
(q/p) are related to, the leading term in this expansion is represented by the so-called 
"spectator effect" contributions. 

Recently, spectator effects have been computed at 0{as) and 0{A.QCD/^h) for the B- 
meson and A;, lifetimes [6] -[9] and for the Cabibbo-favoured transition (6s — > cc) contribut- 
ing to AF5 [10, 11]. In the case of AF^, the complete set of the h-Qco/'mb contributions has 

^Note that the definition of AF^ used in this paper has an opposite sign with respect to one used in 

the BaBar paper. In addition, the measurement itself has a sign ambiguity that can be removed only by 
measuring the sign of cos 2/3, where (3 is the phase of the mixing. This sign is not measured yet but it is 
known to be positive in the SM, see e. g. [3] . 
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been presented in [12], but the QCD next-to-leading order (NLO) corrections have been 
only estimated in that paper in the limit of vanishing charm quark mass, using the results 
of [11]. 

The main result of this paper is the calculation of NLO QCD corrections to in- 
cluding the contributions from a non-vanishing charm quark mass. Our results for the 
NLO Wilson coefficients of agree with ref. [11] and in addition we have computed the 
Cabibbo-suppressed contributions. We have also checked the computation of the Kqcnl^b 
corrections to Tf^^ and T^i at the LO in QCD and the results agree with refs. [10, 12]. Fi- 
nally, we have computed the matching between the QCD and HQET AS = 2 operators 
at the NLO in QCD and found agreement with the results presented in ref. [13]. As phe- 
nomenological applications of our calculations, by using the lattice determinations of the 
relevant hadronic matrix elements [13], we estimate the width differences AF,^, AF^ and 
the CP violation parameters {q/p)d and {q/p)s- 

The inclusion of NLO corrections is important in order to match the scale and scheme 
dependence of the renormalized operators entering the OPE with that of the Wilson co- 
efficients, thus reducing the theoretical uncertainties. In addition, the NLO corrections 
are found to be generally large. Taking into account the finite value of the charm quark 
mass is also an important issue for a proper NLO estimate of AF^ and of the CP viola- 
tion parameters \{q/p)d\ and \{q/p)s\- In the limit rric — > 0, AF^ becomes independent of 
the phase of the mixing amplitude (2/3 in the SM) and CP is not violated in the mixing, 
i.e. \{q/p)d\ and \{q/p)s\ are simply equal to one [14]. Therefore, the dependence on the 
mixing phase as well as the CP violating terms disappear from the NLO corrections in the 
vanishing charm mass limit. 

Our theoretical estimates of the width differences, obtained by using the lattice deter- 
minations of the relevant hadronic matrix elements [13], are 

AF,/F, = (7.4 ± 2.4) x lO^^ , AFrf/F^ = (2.42 ± 0.59) x 10"^ . (4) 

These results can only be compared at present with the experimental determinations given 
in eqs. (1) and (2). The predictions are consistent with the measured values but, given the 
large experimental uncertainties, it is certainly premature to draw any conclusion. Our 
estimates are accurate at the NLO in for the leading term in the l/m^ expansion and 
at the LO for the first power-correction. 

It is interesting to consider also the estimate of the ratio AF^/AFs, for which we obtain 
the prediction 

AFd/AF, = (3.2 ± 0.8) x 10"^ . (5) 

In this ratio the uncertainties coming from higher orders of QCD and AqcD/iTib corrections, 
as well as from the non perturbative estimates of the hadronic matrix elements, cancel out 
to some extent. 

Finally, our predictions for the CP violation parameters are 

[(g/p)^[ - 1 = (2.96 ±0.67) X 10-\ |(g/p)J - 1 = -(1.28 ± 0.27) x 10"^ (6) 

The value of \{q/p)fi\ can be compared with the preliminary experimental determination 
given in eq. (3). Improved measurements are certainly needed to make this comparison 
more significant. 
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We conclude this section by presenting the plan of this paper. In sect. 2 we review 
the basic formalism of the Bq — Bg mixing, introducing all the physical quantities we are 
interested in. The details of the NLO calculation are presented in sect. 3, along with the 
calculation of the AQcn/iTib correction. In sect. 4 we present the theoretical predictions for 
the observables under consideration in both the Bd — B^ and Bg — Bg systems. Finally, the 
analytical expressions for the matching coefficient functions are given in appendix A, while 
detailed definitions of the renormalization schemes used in the calculation and results for 
the matching between QCD and HQET operators can be found in appendix B. 



2 Bq — Bq Basic Formalism 

The neutral B^ and Bg mesons mix with their antiparticles leading to oscillations between 
the mass eigenstates. The time evolution of the neutral mesons doublet is described by a 
Schroedinger equation with an effective 2x2 Hamiltonian 

■ d_ (Bq\ 

'dt \bJ 

The mass difference Am^ and the width difference AFg are defined as 

Am, = m|,-mi, AF, = Fi - F|, , (8) 

where H and L denote the Hamiltonian eigenstates with the heaviest and lightest mass 
eigenvalue respectively. These states can be written as 

The phase of {q/p)q depends on the phase convention of the B states and hence it is not 
measurable by itself. In this paper, we are interested in \{q/p)q\ only. 

Theoretically, the experimental observables Am^, ATg and \{q/p)q\ are related to Mfi 
and F21 in eq. (7) as follows ^ 



Ml, 



"J 
11 

i 91 



"J • 
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Bq 
Br. 



(7) 



(Am,)^--(AF,)2 = 4|M|,P-|F1,P, 



Am, AF, = -4 Re (M|; Tl,) , (g/p), 



2M|i - iV'' 



21 



2M|; - iV% 



21 



(10) 



In the Bd — B^ and Bg — Bg systems, the ratio Fgi/Mli is of 0{ml/m^) ~ 10 ^. Therefore, 
by neglecting terms of 0{ml/m'^), eq. (10) can be simply written as 



Arriq = 2 



211, AF, = -2|M|,|Re 



F-^ 



21 



Ml, 



1 + - Im 
2 



F^ 



21 



Mil. 



(11) 



^For more information about the basic definitions in the 84— B^- and B^— B^-mixing, see for instance [1] 
and references therein. 
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The matrix elements Mfi and Fji are related to the dispersive and the absorptive parts 
of the AB = 2 transitions respectively. In the SM, these transitions are the result of 
second-order charged weak interactions involving the well-known box diagrams. 

The quantity Mli has been computed, at the NLO in QCD, in ref. [15] and it is given 

by 



(12) 



where Xt = rh^/M^, rjB is the QCD correction factor and 5*0 is the Inami-Lim function. 
Here and in the following we use the notation {qq)v±A = 97/^(1 =t 75)5 and {qq)s±p — 
q{l ± 75)5. A sum over repeated colour indices is always understood. 
The matrix element F^i can be written as 

^ ^Di^^^^^l^ / ^ {KfT'i^)KfT\^)) W , (13) 



where "Disc" picks up the discontinuities across the physical cut in the time-ordered prod- 
uct of the effective Hamiltonians. 

The AB — 1 effective Hamiltonian relevant to 6 — > d and 6 — > s transitions is 



n 



AB=1 
eff 



Gp 

71' 



KtKd (CiQi + C2Q2) + v:,Vcd {CiQl + C^Ql) + {c^u) 



u=3 



+ 



d 



+ h.c. 



(14) 



The Ci arc the AB = 1 Wilson coefficients, known at the NLO in perturbation theory [16]- 
[18], and the operators Qi are defined as 



Ql = {biCj)v-A{Ujdi)v-A, 

Ql = {biCj)v-A{cjdi)v-A , 
Q3 = {bidi)v-A^{qjqj)v-A , 

Q5 = {bidi)v-A^{qjqj)v+A , 

QsG 



Q2 = {biCi)v-A{Ujdj)v-A , 
Q2 = (biCi)v-A{Cjdj)v-A , 

Q4 = {bidj)v-A^{qjqi)v-A ■ 
g 

Qe = {bidj)v-A^{qjqi)v+A : 



(15) 



^mAa^'''{l-j'>)t^jdjGl, 



87r2 



Due to the large mass of the b quark, the time-ordered product in eq. (13) can be expanded 
in a sum of local operators of increasing dimension [4, 5]. Up to order l/m^, this expansion 
reads 



i 9 



21 



'247rMB„ 
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(16) 



Figure 1: Feynman diagram contributing to T^i at the LO in QCD. 

where 

Ci = {hiqi)v-A{bjqj)v-A , C| = (biqi)s-p{bjqj)s-p , (17) 

and Sf^j^ denotes the sub-leading l/m^ corrections. Thus, the OPE of at the leading 
order is expressed in terms of only two local four-fermion operators, Of and Ol- 

In terms of the CKM matrix elements V*^Vcq and V^lVtq (the latter regulates the M|i 
contribution), the Wilson coefficients cf and the S^^^ in eq. (16) can be written as 

c! = {v^ivt^r or + 2 v:,v,qvivtq {or - dd + (kikj^ {dt + - 2 dd , 
5i,m - {v,m,r sr;^ + 2 v^v^viv,, - <5,7^) + {vxv^f {5-1 + - 2 <5,7^) 

(18) 

where the labels cc, cm and uu denote the intermediate quark pair contributing to r2i, 
see the diagram in fig. 1. The unitarity of the CKM matrix has been used to drop out 

In this paper, we compute the QCD NLO corrections to Df^ , Df^ , taking into 
account the finite value of the charm quark mass. Our results are collected in appendix 
A. The results for are in agreement with a previous calculation [11]. The functions 
jjcu ^ jjuu j^g^yg been previously estimated [12] taking the limit of a massless charm quark 
of Df'^. As stressed in the introduction, this approximation is not fully satisfactory, in 
particular if one wants to investigate the dependence on the phase of the mixing in Ar^/rc; 
and the \q/p\ parameters. 

As far as the 1/m contributions arc concerned, ^^m, ^^/m' KJm i^^) have been 

calculated in [10] and [12] respectively. We have repeated the calculation by expanding in 
l/rrih the relevant Feynman diagrams and found agreement with the previous results. 

3 Calculation of the Wilson coefficients for Ffi 

In this section, we discuss the basic ingredients of the calculation of the NLO QCD correc- 
tions to Fgi, at the leading power in l/m^, and of the calculation of the l/m^ corrections 
at the LO in QCD. We refer the interested reader to refs. [6, 8] for further details on the 
NLO calculation. 
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Figure 2: Feynman diagrams which contribute at the NLO to the matrix element of the 
transition operator T, both in the full theory (Di) and in the effective theory (Ei). The 
diagrams obtained from Ei, E2, Di, D2, D^, Dq, Dt, Dg, Diq and D12 by a 180° rotation 
(keeping fixed the flavour of the internal lines), as well as those containing the self-energy 
one-loop corrections to the external fields are not shown. 

3.1 NLO QCD corrections at the leading order in l/mt 

Although at the leading order in 1/m;, the heavy quark expansion of the width difference 
could be expressed directly in terms of QCD operators, we found it easier to perform first 
the expansion in terms of HQET operators and then translating the results into QCD, once 
IR divergences are canceled out. 

In order to compute the Wilson coefficients of the AB = 2 operators at the NLO, 
we have evaluated in QCD the imaginary part of the diagrams shown in fig. 2 {full 
theory) and in the HQET the diagrams E^ {effective theory). The external quark states 
have been taken on-shell and all quark masses, except rUb and rUc, have been neglected. 
More specifically, we have chosen the heavy quark momenta pi = ml in QCD and kh = 
in the HQET while, for the external light quarks, Pq = in both cases. We have performed 
the calculation in a generic covariant gauge, in order to check the gauge independence of 
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the final results. Two-loop integrals have been reduced to a set of independent master 
integrals by using the TARCER package [19]. 

IR and UV divergences have been regularized by D-dimensional regularization with 
anticommuting 75 (NDR). As discussed in detail in ref. [6], in the presence of dimensionally- 
regularized IR divergences the matching must be consistently performed in D dimensions, 
including the contribution of (renormalized) evanescent operators [20]. All the evanescent 
operators entering the matching procedure are collected in appendix B. 

As a check of the perturbative calculation, we have verified that our results for the 
Wilson coefficients satisfy the following requirements: 

• gauge invariance: the coefficient functions in the MS scheme are explicitly gauge- 
invariant. The same is true for the full and the effective amplitudes separately; 

• renormalization-scale dependence: the coefficient functions have the correct logarith- 
mic scale dependence as predicted by the LO anomalous dimensions of the AB — 2 

and AB = 1 operators; 

• renormalization- scheme dependence: we performed the calculation in two different 
NDR-MS schemes for the AB = 2 operators (see appendix B and refs. [11, 21] for 
a detailed definition of these schemes) and verified that the NLO Wilson coefficients 
obtained in the two cases are related by the appropriate matching matrix; 

• IR divergences: the coefficient functions are infrared finite; 

• IR regularization: in the limit rUc — > 0, we have also checked that a calculation using 
the gluon mass as IR regulator gives the same results. In this case, the matrix elements 
of the renormalized evanescent operators vanish and hence these operators do not give 
contribution in the matching procedure; 

• comparison with the results in [11]: our results for Df^ defined in eq. (18) and given 
in appendix A agree with those obtained in ref. [11]. 

The AB = 2 effective theory is derived from the double insertion of the AB = 1 
effective Hamiltonian. Therefore, the coefficient functions D'^^, D'^ and D^^ of the AB = 2 
effective theory depend quadratically on the coefficient functions Cj of the AB — 1 effective 
Hamiltonian. The dependence on the renormalization scheme and on the scale of the 
AB — 1 operators actually cancels order by order in perturbation theory against the 
corresponding dependence of the AB = 1 Wilson coefficients Cj. Therefore, the coefficient 
functions only depend on the renormalization scheme and on the scale of the AB = 2 
operators. 

The analytical expressions of the Wilson coefficients are collected in appendix A. For 
illustrative purposes, the corresponding numerical values, at the reference scales Hi = fi2 — 
ruh are presented in table 1. In this table, the LO and NLO coefficients are shown together 
with the pure 0{as) contribution, both in the case of a finite charm quark mass and in 
the charm massless limit ^. The NLO coefficients refer to the QCD operators renormalized 

^Note that, in what we call the pure 0{as) contribution in table 1, the NLO corrections to the AB = 1 
Wilson coefficients are not included. For this reason, the NLO results in the table 1 differ from the sums 
of the LO coefficients and the pure 0{as) contributions. 
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LO 


0{a,) 


NLO 


e»(«,) 


NLO 


















0.260 


-0.123 


0.174 


-0.123 


0.174 


JJUU _ J^CU 




0.103 





0.105 


-0.034 


0.072 


Df + Df - 




0.004 





0.004 


-0.005 


-0.001 


r)uu 




-1.56 


0.44 


-1.00 


0.44 


-1.00 






-0.035 





-0.033 


0.012 


-0.022 


+ D^" - 


2L)f 


0.016 





0.015 


-0.014 


0.001 



Table 1: Values of the Wilson coefficients defined in eq. (18) computed at the reference 
scales Hi = 1^2 = rub = 4.75 GeV. The LO (first column) and NLO (third and fifth 
column) coefficients are shown together with the pure 0{as) contribution (second and fourth 
column). The results given in the second and third column are those obtained by neglecting 
the 0{asmc) corrections. 

in the NDR-MS scheme of ref. [11]. In order to obtain the numerical results we used the 
central values of the input parameters given in table 2. 

By looking at the results shown in table 1, we see that the NLO corrections are large 
and the effect of charm contributions is crucial to estimate accurately all the terms in 
eq. (18). 

3.2 l/rrib Corrections to Ffi at the LO in QCD 

In order to compute the l/mj, correction at the LO in QCD, the imaginary part of the 
diagrams in fig. 1 have been evaluated between on-shell quark states and expanded in the 
quark momenta ^. The result of the full theory has been then matched at 0{l/mb) onto 
the following independent set of QCD operators, 

'rri — — 1 — i — 

Rl = A(6iSp(l - 'y,)D%){b,{l - 75)g,) , Rl = —{k{l - -i,)iD,qi){^Y{'^ - ^,)q,) . 
nijj rrib 

(19) 

"'In the heavy quark expansion we consider the strange quark mass of 0{Aqcd), whereas m„ and 
are neglected. 
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Since derivatives acting on b fields scale as mj, the above operators are manifestly of order 
l/rrih. The results for the l/m^ terms defined in eq. (18) read 

= VT^z{{l + 2z)[K,{{Rl) + 2{Rl))-2K,{{Rl) + {in))] 

- ^^jK,{{Rl)+2{Rt)) + 2K,{Rl)]), 

^Til = (l-^)^((l + 2z)[i^2(W + 2(i?l))-2i^a(^?) + (^l))] (20) 

- ^^(i^i((i?2^) + 2(i?I)) + 2ir2(i?I))), 
<^V^ = [i^2(W + 2(i?I))-2i^i((i??) + (i?l))], 

where Kx and K2 are related to the AS = 1 coefficients, = 3 +2 Ci C2 and ^2 = C*!) 
and 2; = {mc/mijY. Our results agree with previous calculations [10, 12] ^. 

4 Width differences and CP violation parameters in 
Bq — Bq mixing 

In this section we present the theoretical estimates for the quantities of interest in this 
paper: the width differences AFd/Fd and AFs/Fs and the CP violation parameters \ {q/p)d\ 
and \{q/p)s\- 

Our estimates neglect terms of 0{al), 0{as {kqcD/^b)) and 0{{KQCD/^bY)-, whereas 
the charm quark mass contributions have been fully taken into account. For the relevant 
hadronic matrix elements we have used the lattice determination of ref. [13]. 

In eq. (11), the width difference AFg and the parameter \{q/p)q\ for neutral B mesons 
are expressed in terms of the real and the imaginary part of F^i/Mli respectively. Taking 
into account the expressions obtained in eq. (12) for Mfi and eqs. (16)-(18) for F21, we can 
write 



AF„ 



A///y /C 



X 



{Bq\Ol\B,) 



,i=l,2 



^uu , cos2^ 



+ 



cuuq I 



cos2 13 q 

Rtq 



(Df + D]^"-2Dn -2 

cos /3q 



COS (3q 



R 



uu 
i 



■tq 



(cccq I sruuq o fCU(j\ o COa fjq / ruuq rctigX 

^tq 



(21) 



q 

.p. 



- 1 



k: 



2{Bq\Ot\Bq) 



X 



^For the reader's benefit, we observe that the expression ofSy^ in eq. (49) of ref. [12] has been expanded 
up to 0{z^). 
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E {B,\Ot\B,) 

\i=l,2 



sin 2(3q , ^ „ T^cu\ o /^g 



(Df + D]'" - 2 ) - 2 



+ 



sin 2/?y 



/m 



sin/?g 
R 



tq 



ruu q ecu q 



(22) 



The two operators entering these expansions at the leading order, 0\ and Ol, are defined 
in eq. (17). The coefficient K, is given hy K, — 4:7rml/{3My^rrjBSo{xt)). In addition we have 
used 

V:,Vj{V,lVt,) = -e^"^/Rt,. (23) 

f3d and Rfd are the usual angle f3 and the side Rf of the unitarity triangle respectively, 
whereas Ps and Rts parameterize the Cabibbo-suppressed contributions to the Bg system. 
For completeness, we give their expansion in terms of the parameters of CKM matrix, up 
to and including C(A^) terms. They read 



Rtd = pT + V^ , sinPd ^v/Rtd, 

Rts^l + + iA^(2p + Tj^) , sin = -A^ry [1 + A^ (1 - p)] . 



(24) 



The Cabibbo-suppressed contributions are practically irrelevant for AF^ so that, to an 
excellent approximation, one can put Ps — and = 1 in eq. (21) to obtain 



(25) 



The contributions neglected in the previous equation, however, give rise to the CP violation 
effects in the Bg — Bg mixing, which are accounted for by the deviation of the parameter 
\{q/p)s\ from unity. Using f3s = and i?^^ = 1 in eq. (22) gives in fact |((//p)s| = 1. 

For the matrix elements entering our calculation (see eqs. (17) and (19)), we use the 
following parameterization 



{B,\Ot\B,) = -flMlBl 
(BMW^l'^JlMlBl^, 



{B,\Ol\B,) = 



5 PeM 



B'. 



3 (rrib + niq^'^ ^ ' 



Ml 



{B,\Rl\B,)^--flMl 
l) 5^3, {B,\Rl\B,)^-fim 



'Ml 



mt 



'Ml 



mt 



{B,\Rl\B,) = -flMl 

(26) 

Among these 5-parameters, Bf and i^f are the most widely studied and well known 
in lattice QCD [13], [23]- [29].^ In this paper we use the results of ref. [13], in which 
the complete set of AB = 2, dimension-six, four-fermion operators has been determined 
in the quenched approximation of QCD. For B^ and B2, the results of [13] are in very 



^For estimates of these matrix elements based on QCD sum rules, see refs. [30]- [32]. 
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good agreement with those obtained in [24, 25] by using the lattice NRQCD approach. In 
addition, it has been shown in refs. [28, 29] that the effect of the quenching approximation 
for these quantities is practically irrelevant. Thus, the systematic uncertainties in the 
present lattice estimates of Bf and are quite under control. 

To our knowledge, the matrix elements of the operators Rf defined in eq. (26), entering 
the l/iTLii corrections, have been only estimated so far in the VSA. Using the complete 
set of operator matrix elements calculated in [13], however, two of the four independent 
parameters S^. can be also evaluated. Besides the operators Of 2: the complete basis 
studied in [13] includes 

^3 = (biqj)s-p(bjqi)s-p , Ol = {biq.i)s-p{bjqj)s+p , Ol = (biqj)s-p(bjqi)s+p , (27) 
whose matrix elements are parametrized as 

(B,\Ol\B,) 
(BM\B,) 
{B,\Ol\B,) 

We then notice that the operator R\ defined in eq. (19) is trivially related to 0\, 

Rl^Hhol (29) 
nib 

In addition, by using the Fierz identities and the equations of motion, the operator Rl can 
be expressed as [10] 

2Rl = Ol + Ot/2 + Ol - 2^ Ol + i?f . (30) 

nib 

Thus, eqs. (29) and (30) can be used to get rid of the two parameters i?^^ and 5^^. 

The values of the S-parameters obtained in ref. [13] and used in our calculation are 
collected in table 2. Concerning the unknown matrix elements of the operators R^ and 
i?3, they have been estimated in the VSA and we have included a 30% of relative error to 
account for the corresponding systematic uncertainty ^. 

In order to obtain the theoretical predictions presented in this paper, we have performed 
a Bayesian statistical analysis by implementing a simple Monte Carlo calculation. The 
input parameters have been extracted with fiat distributions, by assuming the central 
values and standard deviations given in table 2. 



4.1 Width differences: ATa/Td, Ar^/r^ and AT a/ AT ^ 

The theoretical predictions for the width differences Ar^/ra and AF^/Fs, as obtained from 
eq. (21) as functions of the CKM matrix elements and of the S-parameters entering at the 

''Since terms proportional to nid/mi, are neglected in our calculation, the values of Bf and are 
presented in the table only for completeness. 



11 



mb = 4.75 ±0.11 GeV 
mt{mt) = 165.0 ±5.0 GeV 
Mw = 80.41 GcV 
t{Bs) = 1.461 ± 0.057 ps 
Mb^ = 5.279 GeV 
r7fe(mb) = 0.85 ± 0.02 


mjmb = 0.30 ± 0.02 
ms{mb) = 87±21 MeV 
Arud = 0.502 ± 0.006 ps~^ 
T{Bd) = 1.540 ±0.014 ps 
Mb, = 5.369 GeV 
i = 1.24 ±0.04 


A = 0.2240 ±0.0036 

/7 = 0.162 ±0.0 16 


a,{mz) = 0.118 
ri = 0.3 17 ± 0.027 


B{ = 0.87 ± 0.05 
Bl = 0.84 ± 0.04 
Bl = 0.91 ± 0.08 
Bl = 1.16 ±0.07 
Bl = 1.75 ±0.20 


Bf/Bf = 0.99 ± 0.02 
B^/B^ = 1.01 ± 0.02 
Bl/Bl = 1.01 ± 0.03 
Bl/Bf = 1.01 ± 0.02 
B^/Bl = 1.01 ± 0.03 



Table 2: Central values and standard deviations of the input parameters used to obtain the 
theoretical estimates of the width differences and of the CP violation parameters. When the 
error is not quoted, the parameter has been kept fixed in the numerical analysis. The values 
of rrib and ni,, refer to the pole masses, while Wis Q^nd fn,t are the masses in the NDR-MS 
scheme. The B-parameters are renormalized in the NDR-MS scheme at the scale rUb- The 
definition of the renormalization scheme can be found in [11] for Bl, and B^ and in 
[20] for Bl and Bl, see also appendix B. 



LO in l/ruf,, can be expressed as 



= 10" 



1 



0.96(10) ±5.48(34)^ -2.9(1.3)^ 

cos/3 / B^ 1 



Bf 



Bf 



Rt 

COS 2/? 



0.80(13) ± 0.24(7)--f - 0.03(6) 

B. 



Bf 



10 



-2 



{XRt 



-0.005(4) - 0.007(16)-| - 0.008(17)-^ 

Bl Bl , 



0.02(2) ±0.77(7)|i- 0.47(18)-^ 



(31) 
(32) 



These formulae are one of the main results of this paper. In the above expressions, the 
terms proportional to l/Bf represent the contributions of the l/rub corrections. Note also 
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Figure 3: Theoretical distributions for the width differences in the Bd and Bg systems, 
predictions are shown at both the LO (light/red) and NLO (dark/blue). 



The 



that, in order to obtain the prediction for AFs/Fs, the mass difference Am^ has been 
evaluated in terms of Arud by using 



Arris = Amd 



Mb, 



V, 



ts 



td 



(33) 



where e=(/B.v^)/(/B,/Bf^ 

The errors on the numerical coefficients presented in eqs. (31) and (32) take into account 
both the residual NNLO dependence on the renormalization scale of the AB = 1 operators 
and the theoretical uncertainties on the various input parameters. To estimate the former, 
the scale fii has been varied in the interval between mb/2 and 2mh. These errors are 
strongly correlated, since they originate from the theoretical uncertainties on the same set 
of input parameters. For this reason, they have not been used to derive our final predictions 
for the width differences. For these predictions, we quote the average and the standard 
deviation of the corresponding probability distribution functions obtained directly from 
the Monte Carlo simulation, namely 



ATd/Td = (2.42 ±0.59) x 10 



-3 



AF,/F, = (7.4 ±2.4) x 10~ 



(34) 



The theoretical distributions are shown in fig. 3. Notice that the difference between 
the NLO and LO distributions is remarkable. The NLO corrections decrease the values of 
both ATd/Td and AF^/F^ by about a factor two with respect to the LO predictions. We 
find that the total error on the width differences is dominated by the uncertainty on the 
6-quark mass in the first place, followed closely by the one due to the renormalization scale 
variation. Other contributions to the error, coming from the ratio rric/mb, the B- and the 
CKM-parameters are smaller, though not entirely negligible. 

In fig. 4 we show the dependence of ATd/Td on the mixing phase (3 as given in eq. (31). 
NLO corrections in the limit of vanishing rric just shift the LO curve, while mc-dependent 
terms modify also its profile. In principle, a measurement of ATd/Td could allow a deter- 
mination of (3, but we find that the dependence is so mild that the extracted value would 
be affected by a very large error. 
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Figure 4: Prediction for AF^/F^ as a function of (3 as obtained at the LO and NLO. The 
NLO estimate o/AF^/F^ at the measured value of (3 is also shown. 



Another interesting prediction concerns the ratio AF^/AFs. This quantity is of partic- 
ular interest, because the uncertainties coming from higher order QCD and AqcD/fnb cor- 
rections, as well as those coming from the non-perturbative estimates of the B-parameters, 
are expected to cancel in this ratio to some extent. From our numerical analysis, we obtain 
the NLO prediction 

AFrf/AF, = (3.2 ± 0.8) x 10^^ _ (35) 

The corresponding theoretical distributions at the LO and NLO are shown in fig. 5. As 
can be seen from the plot, this quantity is practically unaffected by the NLO corrections. 



4.2 CP Violation parameters: \{q/p)d\ and \{q/p)s\ 

CP violation in the mixing shows up only for a non-vanishing value of the charm quark 
mass, as can be explicitly verified from the expression of \{q/p)q\ in eq. (22). In addition, 
it is apparent from eq. (22) that \{q/p)s\ — 1 is suppressed with respect to \{q/p)d\ — 1 by 
a factor sin /5s ~ A^. 

For the parameters \ {q/p)d\ and \ {q/p)s\, expressed in terms of the CKM matrix elements 
and the 5-parameters Bf, B2, we obtain the expression 



10" 



sin P, 



R, 



tq 



^ (5.19(84) + 1.57(47)|| 



0.22(39) 
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Figure 5: Theoretical distribution for the ratio ATa/ATs as obtained at both the LO 
(light/red) and NLO (dark/blue). 




Figure 6: Theoretical distributions for \ {q/p)g\ — l in the and Bg systems. The predictions 
are shown at the LO (light/red) and NLO (dark/blue). 



sin 2l3q 



B'^ 

0.04(3) + 0.04(10)-! + 0.05(11)-^ 



B1 



Bl. 



(36) 



where q = d,s. 

By using the values of the CKM and i?-parameters given in table 2 we obtain from the 
Monte Carlo analysis the final predictions 



1 = (2.96 ±0.67) X 10"^ 



1 = -(1.28 ±0.27) X 10" 



(37) 



The theoretical distributions for these quantities are shown in fig. 6 both at the LO and 
NLO. One can see that the effect of the NLO corrections turns out to be important also 
for these quantities. The error on \{q/p)q\ — 1 is largely dominated by the uncertainties 
on rric/mh and on the CKM parameters. We note in particular that the dependence on 
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Figure 7: Prediction for the CP violation parameter \{q/p)d\ — 1 as a function of j3 as 
obtained at both the LO and NLO. The NLO estimate of \{qlp)d\ — I at the measured value 
of P is also shown. 

the 5-parameters is practically negligible, making these predictions free from hadronic 
uncertainties. 

Finally, we show in fig. 7 the LO and NLO predictions for \{q/p)d\ — 1 as functions 
of p. In this case, NLO corrections without charm mass leave the LO result practically 
unchanged, while the mc-dependent terms give a sizable contribution. The dependence on 
f3 is quite strong, allowing for a determination of the mixing angle once \{q/p)d\ will be 
measured. 



5 Conclusions 

The main result of this paper is the NLO QCD calculation of the Wilson coefficients 
entering the heavy quark expansion of V^i and r2i, including the effects of a non- vanishing 
charm quark mass. For we find agreement with a previous result [11]. Using our 
results, we improve the theoretical predictions for different observables in the neutral B- 
meson systems, namely the width differences ATd/Td and AFs/F^ and the CP violation 
parameters \{q/p)d\ and \{q/p)s\- 

Using formulae including the NLO corrections at the lowest order in the l/m^ ex- 
pansion and the first power corrections at the LO in as, we find 

AF AF AF 

— ^ = (2.42 ±0.59) X 10-^ — ^ = (7.4±2.4) X 10-^ — ^ = (3.2 ± 0.8) x 10"^ 
r<i Ts AF^ 
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- 1 = (2.96±0.67) X 10-^ |(g/p)J - 1 = -(1.28 ± 0.27) x 10-^ (38) 

NLO QCD corrections are in general important for theoretical consistency and give 
sizable contributions to the Wilson coefficients considered in this paper. We find that the 
charm quark mass effects at NLO are numerically important, in particular for \{q/p)q\-, and 
should be included in the phenomenological analyses. 
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Note added 

Just before this paper was finalized, a paper by Beneke et al. on the same subject was 
submitted to the e-print archive [34]. Comparing the results, we found that the NLO 
contribution to the Wilson coefficients D"^ of eqs. (18) and (39) differ from the ones given 
in eq. (25) of ref. [34]. Had we assumed that the diagrams Dq, D7 and Ds in fig. 2 were 
equal to the corresponding "rotated" ones {i.e. those diagrams where the gluon is attached 
to the other external quark leg with the same flavour and to the other virtual quark line), 
we would have obtained the result of ref. [34]. However this assumption, which is correct 
when the two virtual lines correspond to the same flavour, docs not hold when the two 
quarks are different. This may explain the origin of the discrepancy. 



Appendix A: Analytical Results for the AB = 2 Wilson 
Coefficients 

In this appendix we collect the analytical expressions of the Wilson coefficients, both at 
the LO and NLO. The LO coefficients have been computed in refs. [10, 12] and are given 

here for completeness. 

The coefficient functions -D^", -D™ and D'^'^ defined in eq.(18) depend quadratically on 
the coefficient functions Cj of the AB = 1 effective Hamiltonian and can be written as 

Drif^2) - E ^^.(/^l)^.(/^l)^M.(/^l,/^2) + /C2(/.l)2p,^^2(/^l,/^2) 

i=l,2r=3,6 

+ TZ ^2 CsG (-Pfe,28 + Pk,2s) + E QCr {Pk,ir + Pk,ir) ' 

i=l,2r=3,6 
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mi^2) = E a(//l)C,(/.i)F,%(/^l,/^2) + ^C2(/.l)2p,%(//l,//2) 

+ 2-^C2CsGPk,2S + 2 E EC^CrP^,^r■ (39) 

i=l,2r=3,6 

The functions F^\j arc obtained from the insertion of the operators Qi and Qj in the 
Feynman diagrams Di — Diq of fig. 2, whereas the diagrams Du and contribute to the 
functions -P^ 22(a*1' A*2) and P^^^ respectively. Contributions with the double insertion of 
penguin operators have been neglected, since the Wilson coefficients C^-Cq are numerically 
small. 

We distinguish the leading and next-to-leading contributions in the coefficients F^'^-j by 
writing the expansion 

+ (40) 

where {qq') = {{uu), (cu), (cc)}. Since only the sum of P^i2 and contributes to 

eq. (39), we just give in the following the average of the components 12 and 21. Moreover, 
we do not write exphcitly the results for the functions since they can be obtained by 
taking the limit mc —> of D^. 

In terms of the ratio z — ml /ml, the LO coefficients A^^^- read 



^Mi = i (2-3^ + ^3^ , =^Jr2i = 5 (2-3^ + ^3) , Al^^^ = \{l-zf , 

1 + 

(41) 



^241 = 3 (1 - 0)' (1 + 2 ^) , Ai^i2 = ^^21 = (1 - (1 + 2 ^) , Ai^22 = -(1 - ^)' (1 + 2 z) 



Ai^ii =3^1^^(1 + 2^) , ^§^12 = (1 + 2 0) , ^i^22 = -vT^4^(l + 2 0). 

(42) 

The NLO results for the coefficients Bf^^- are presented in the NDR-MS scheme of 
ref. [21] for the AB = 1 operators and the NDR-MS scheme of ref. [11] for the AB = 2 
operators, in QCD. We find 

^Mi = ^-37z+^ + ^ + 2(l-z)2 (5 + z)logX2-4(l-z)^ (5 + 7z)log(l-z)- 

' D 2 6 

2 z (10 + 14 z - 15 z2) log z + 8 (2 - 3 z + z^) log(l - z) log z + 16 (2 - 3 z + z^) Li2(z) , 
-B^ii = (10-33 z + 54z2_ 31^3) -8(1- z)^ (4+ 14z-3z2) log(i_^) + 
8 z (2 - 23 z + 21 z^ - 3 z^) log z - 

16(l-z)' (l + 2z) (2 1oga;2-log(l-z) logz - 2Li2(z)) , (43) 

-Bfri2 + -S?r2i -502 + 912z-387z^-23z3 ^■^n7^A^^ ,2,,^,, 

— ■ — ^ — = ^ (1-^) (17-|-4z) loga;i-|--(l-z) (5 -|- z) log a;2 - 

^^^^1^ (2 + 33z + 94z2) log(l-z)-^(80 + 69z-126z2) logz + 
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-{2-3z + z^) (log(l-^) log^ + 2Li2(-2)), 

^ii^— — ^ = —- ^— ^ (1 + 15^ + 470^-12^3^ log(l-^) + 

2 9 o z 

2 

-z (8-93z + 87z2- 12 2^) logz- 

(1 + 2 0) (3 1oga;i+4 1ogX2-2 1og(l-0) log^ - 4^2(0)) , (44) 

r . -136- 159 + 738 02-443 03 
^^22 = -y (1- 50 + 402) + 2(1-0)' (5 + 40) logxi + 

-{l~zf (4-0) loga;2 + ^-^^ (7 + 32 0^ + 303) log(l - 0) - 

- (62 - 390 -3002 + 3^3) log0+ ^ " ^ - + ^ (log(l - 0) log0 + 2Li2(0)) , 
6 3 

8 28 
^2^22 = ^7r2 (1 + 0-202) -— (5 + 30- 2702 + 1903) - 16(1 - 0)^ (1 + 20) loga;i + 

— (1-0)' (I + 20) \ogX2-^^^—^ (1 - 120- 16 02 - 303) log(l-0) + 

3 30 ^ ' 

\z (2 - 30 + 18 02 - 303) log0+^ (1 - 30- 602 + 803) (log(l - 0) log0 + 2 ^2(0)) j(45) 

VI - 4 (109 - 226 + 168 02) , „ 
Bl\, = !^ ■ ^ - (52- 1040- 1602 + 5603) logcr + 



2 (5-80) Vl-40log.T2-12Vl-40 (3-20) log(l-40) + 4 (I3-IO0) Vl - 4 log - 
16 (1 - 30 + 202) (-siogS^ + 2 logcr log(l - 40) - 3 log cr log0 + 4Li2(cr) + 2Li2(cr2)) , 
-8 VI -40 (5- 230-420 



2\ 



= ^ i-- 16 (4 - 20 - 702 + 1403) logcr -32V1- 40 (I + 20) loga;2 



48 Vl - 40 (1 + 20) log(l - 4 0) + 64 Vl - 40 (1 + 2 0) log + 

16 (1 - 4 02) (31og2a + 2 1ogalog(l-4 0)-3 1ogalog0 + 4Li2(cr)+2Li2(cr2)) , (46) 



_ - VI - 4 (127 - 199 - 75 02) (2 - 259 + 662 02 - 76 03 - 200 0^) log ct 
^^•12 - g + ■ YTz 

(17-260) yT34llogxi + ^^^^M^JHli!^_4Vr^ (3-20) log(l-40)- 

VI -40 (2 - 255 + 316 02) log 

120 ^ 

16 (1 - 3 + 2 02) (3 log2o- + 2 log (7 log(l - 4 0) - 3 log (7 log + 4 Li2(cr) + 2 Li2(cr2)) 

-2V1 - 40 (68 + 49 - 150 02) 2 (l - 35 + 402 + 7603 - 1000^^) ^^g^ 
^2.12 = g + + 



(16 - 64 02) log2o- - 8 VI -40 (1 + 2 0) log xi 



32 VI- 4: z (1 + 20) log 2:2 



2 VI -40 (1 - 33 0- 76 02) log0 



16V1 -40 (1 + 20) log(l-40) !^ ^ — — + 

3 

16 (1 - 4 02) (2 log (7 log(l - 4 0) - 3 log C7 log + 4 Li2(<7) + 2 02(0-2)) 

3 ' 
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(47) 



-7r2(l-10z) VI -4z (115 + 632 z + 96^2) (7 + 13z - 194^^ + 304^3 _ 64^-*) loga 



22 



18 6^ 



2^/n^(5-2 0)loga;i + ^^^ ^"^"^ ^ - 4 (1 - 6 z) VT^4^ log(l - 4^) + 

o 

(l3-54z + 8z2) logo- log(l-4z) y/l-Az (7 + 27^- 250^2) bg^ 
3 6^ 

^7 -^9 ^4 2w, 2 , , , , 4 (5-12^ + 4^2) ^ 4 (4-21.3 + 2.^2) Li2(a2) 
(7 - 32 z + 4 z j (log o- - log (J log z) H 1 

O O 

STT^il + ^z) 8V1 -4z (l9 + 53z + 24z2) 4 (l + 7z + 10 z^ _ 68 z^ + 32^4) log cr 

2,22 — o ' 



3z 

8 (1 + 2 zf lo^a - 16 VT^z (1 + 2 z) log + 32 x/T^ (1 + 2 z) l og^ ^ 

o 

, , . , . 8 (l + 6z + 8z2) loga log(l-4z) 

16\/l-4z (l + 2z) log(l-4z) !^ ^—5 i 

4^/^^ (l + 9z + 26z2) logz , ^ , ^ , ^ 
|-8(l + 2z) logo- log z + 

3 z 

32 (l-4z2) Li2(a) 32 (l + 3z + 2z2) Li2(a2) 



3 

where a is the ratio 



(48) 



a = ^-^"^ (49) 

and we have defined xi — \x\lm\, and X2 = 1^2/ i^b- 

The contributions of the diagram Du in fig. 2 and of the insertions of the penguin and 
chromomagnetic operators, read 



pen 



22 



1 2z loga;i Vl - 4z (1 + 2 z) (2 + 3 logtr + 6 logxi) logz 
'27 ~ "9 9 54 ^ ~W 



^ 8 16z 8 1og.xi 4vT^^(l + 2z) (2 + 3 1oga + 6 1oga;i) 4 logz 

2,22 27 9 9 27 9 ' ^ 

_ -2V1-4Z (l + 8z + 12z2) _ logo- 4z2logc7 16 zMog cr _ 

^'22 ~ 27 9~ ^ 3 ^ 9 
VI -4z (l + 2z) (2 logxi - logz) 

_ 16V1-4Z (l + 8z+12z2) 8 logg 32 z^ logo- _ 128 z^ loga 

2'22 ~ 27 ^ 9 3 9 ^ 



8V1 -4z (l + 2z) (2 logxi - logz) 



Pf,13 = 3 {I-Z) , Pf,23 = VT^ (1 - Z) , Pf,i4 = VT^ (1 - Z) , 

P2V3 = 3\/r^^(i + 2z) , p|,23 = \/r^^(i + 2^) , P2% = vr^^(i + 2z) 

= \ ( 1-4^) ^ , Pf,i5 = 9 z VT^4i , Pf_25 = 3 ^ ^r^^z , 

Phi = -VT^z (1 + 2 z) , P2-,i5 = , P^^^ = . 

pfi6 = 3zvr^^, pf,26 = 3^%/r^^, A% = -^ yr^4 i (i + 2z) , 

^2.16=0, ^1,26 = 0, P|,28 = |\/r^4l(l + 2z). 

The functions P^j^ can be obtained from P^^^- by taking the hmit 2; — > 0. 



(51) 

(52) 

(53) 
(54) 
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Appendix B: Evanescent Operators and Renormaliza- 
tion Schemes 



In this appendix, we define the NDR-MS scheme chosen to renormahze the QCD, AS = 2 
operators whose Wilson coefficients are given in appendix A. This scheme is introduced in 
ref. [11] by giving three prescriptions (cqs. (13)-(15) of that paper) which define impficitly 
the relevant evanescent operators. In order to make it more explicit, we write out in this 
appendix the complete basis of operators and evanescent operators defining this scheme. 
Although this renormalization scheme is the only one needed to define our results, we think 
it may be useful, for future applications, to present also the matching matrix which relates 
the QCD operators in the NDR-MS of [11] to a corresponding set of HQET operators. 
This matrix has been computed and used in the intermediate steps of our calculation. To 
this purpose, wc consider a NDR-MS scheme for the HQET operators which is a simple 
generalization of the one defined in ref. [22]. In addition, since a second NDR-MS renor- 
mahzation scheme has been introduced in ref. [21] for the AS = 2 operators in QCD, we 
also discuss this scheme and present the corresponding NLO matching matrix relating the 
QCD operators to the HQET ones. 

To start with, we introduce the following set of AS = 2 four-fermion operators. 

Of = (biqi)v-A(bj1j)v-A , 01 = {hiqi)s-p{hjqj)s-p , 

(55) 

0| = {hi a'^ qi){bj cr^^Lqj) , 0| = {bi qj)s-p{bj qi)s-p , 

where, for any string F of Dirac matrices, we define Tl = r(l — 75). 

Concerning the evanescent operators, it is worth to recall that they are renormalized, in 
any given renormalization scheme, in such a way that their matrix elements vanish on IR 
finite, physical external states. Consequently, the anomalous dimension matrix elements 
mixing the physical and the evanescent operators vanish to all orders [33]. Furthermore, 
while the evanescent operators usually only enter the definition of the renormalization 
scheme, in some cases they can also contribute beyond the LO to the matching of the 
physical operators. In particular, when the matching is performed in the presence of IR 
divergences, as in our calculation, one should properly take into account the contribution of 
the matrix elements of the evanescent operators to the matching conditions of the physical 
operators [6, 20]. 

We now proceed by defining in details the several renormalization schemes discussed 
above. 



NDR-MS scheme for QCD operators of ref. [11] 

This scheme concerns AS = 2 operators in QCD. The four-dimensional basis involves the 
operators Of, 0\ and 0\ in eq. (55), whereas the evanescent operators are 

SI = {h qj)v-A{bj qi)v-A - 01 , 

SI = YYYl Qj)ibj l.lulpL g.) - (16 - 4 £) {SI + 01) , 
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SI = (6.<g,)(6,a^.Lg.)-(4-2£)Ol-(8-8£)0|, (56) 

SI = (6,7V7^g.)(&,7M7.7pLg,)-(16-4e)0?, 

SI - (6,<g,)(^.-^M-Lg,)-(4-2£)0|-(8-8£)0|. 

NDR-MS scheme for QCD operators of ref. [21] 

This scheme is defined by choosing 01 . O2 and O3 in eq. (55) as operators of the physical 
basis and the following set of evanescent operators 

Ml = {hi qj)v^A{bj qdv-A - Of , 

Ml = (6.arg,)(6,a^,Lg,)- 6 02^ -0,^/2, 

M! = ikq,)s-pibj1^)s-P + 01/2 -01/8, (57) 

Ml = (6,7^7V7Zg,)(&,7M7.7p7.Lg.)-16M|-64M|-640|-16 (l-£) 0|, 

Ml = (6.7V7^g.)(&,7M7.7pLg,) - (16-45) O?, 

Ml = (6,7^7V7£?i)(&i7M7.7p7aLg,)-(64- 96£)0|-(16-8£)Ol. 

NDR-MS scheme for HQET operators 

In the AB = 2 HQET Hamiltonian, we choose the operators Of and Of in eq. (55) as 
operators of the physical basis and the following set of evanescent operators, 

Hi = {hi Qj ) v-Aihj qi)v-A-Ol 

m = {hrYYLqj)ihl,l^lpLq^d-{lQ-^s){H! + Ot) 

HI = {:ha^^q^){h,a^,Lq^) + {^-'2eT){Hl + OD + {A-2ea){Hl-0\/2-Ol) 

HI = 01 + 01/2 + 01 (58) 

HI = {krYYLq^){bJ7,lulpLqJ)-{l6-^e)Ot 

HI = 0,^ + (4-2er)0? + (4-2£(t)0| 

where O3 and 0| are also defined in eq. (55). The parameters a and r are not fixed 
and enter the definition of the renormalization scheme. For the specific choice a — r this 
scheme reduces to the one considered in ref. [22]. In four dimensions the operators Hf — Hl 
vanish because of the Fierz identities, the Chisholm identity 

7^75 (59) 

and the equation 7^6 = b, verified by the 6-field operator in the static limit. Notice that, 
since if 4 and Hi are evanescent, both O3 and O4 in the HQET are reducible in terms of 
Oi and O2 and need not to be included in the basis. 

Matching between QCD and HQET operators 

We now discuss the matching, at the NLO, between QCD and HQET operators. The 
matching condition can be written as 

Cfc(m6)QCD = Cki Oi{mb)HQET + 0{l/mb), (60) 
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where a common renormalization scale fi = mj, has been chosen. Let us stress that the 
equation above is only valid once the operators are sandwiched between external on-shell 
states. 

For the QCD NDR-MS scheme of ref. [11], the matrix C has been given at 0{as) in 
ref. [13], by considering in the HQET the renormalization scheme defined above with the 
particular choice a — t — 0. We have repeated the calculation for generic a and r and 
obtained 



/ Oi(m,) 
02(^5) 
V 04{mb) 



ref. [11] 



QCD 




Oi(mfe) 
02{mb) 



(61) 



HQET 



-14 
39 - 3(7 + 4r 



24 

121 + 3(7 



4r 



-8 

25 -a 

3 

11 + a 



V 



(62) 
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For a = T = this confirms the previous results. 

Notice that in QCD the operator Of is in general independent from the operators Of 
and Of- However, its matrix elements between external on-shell states can be expressed 
in terms of the matrix elements of Of and Of- This relation has been given in ref. [11] in 
terms of the operator 

Rl = 0l + 0t/2 + 0l (63) 

The on-shell matrix elements of Rq are of short-distance origin and their values can be 
easily obtained by using the matching coefficients between QCD and HQET operators 
given in eq. (61). 

When the QCD operators are renormalized in the NDR-MS of ref. [21] the matrix C is 
given by 



/ \ ref. [21] 

Oi(mfe) 



02(mfe) 
Osirrib) 



QCD 



/ 



47r 



Oi(mfe) 
02{mb) 



(64) 



HQET 



35 



-14 
3ct + 4t 



24 

237 + 3a - 4r 



6 



13 



3 

,27 + (J 



(65) 
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By using the above results one also finds that in this scheme the matching of the operator 
O4 onto the HQET operators is given by 



, . xrcf. [21] 



-1/2 -1 ) + 



An 



101 + 3(7 - 4t 7 + a 
24 3 



Oi(mb) 
02{mb) 



/ HQET 

(66) 
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